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Abstract 

A method is suggested for the calculation of the DeWitt-Seeley-Gilkey 
(DWSG) coefficients for the operator \J —V^ + V{x) basing on a generaliza- 
tion of the pseudodifferential operator technique. The lowest DWSG coef- 
ficients for the operator a/— + V{x) are calculated by using the method 
proposed. It is shown that the method admits a generalization to the case 
of operators of the type (— + V^(X))^/™, where m is an arbitrary rational 
number. A more simple method is proposed for the calculation of the DWSG 
coefficients for the case of strictly positive operators under the sign of root. 
By using this method, it is shown that the problem of the calculation of the 
DWSG coefficients for such operators is exactly solvable. Namely, an explicit 
formula expressing the DWSG coefficients for operators with root through 
the DWSG coefficients for operators without root is deduced. 
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1 Introduction 



The algorithms for obtaining the asymptotic heat kernel expansion for second 
order differential operators on a Riemannian manifold are well known [1-3]. 
The most popular is that of DeWitt [1,4] which uses a certain ansatz for 
heat kernel matrix elements. The method possesses the explicit covariance 
with respect to gauge and general-coordinate transformations. However, the 
DeWitt technique does not apply to higher-order operators, nonminimal op- 
erators and, generally speaking, to operators whose leading term is not a 
power of the Laplace operator. Recently, using the Widom generalization 
[5] of the pseudodifFerential operator technique a new algorithm was devel- 
oped [6] for computing the DeWitt-Seeley-Gilkey (DWSG) coefficients. The 
method is explicitly gauge and geometrically covariant and admits to carry 
out calculations of the DWSG coefficients by computer [7]. As was shown 
in [8,9], the method permits a generalization to the case of Riemann-Cartan 
manifolds, i.e., manifolds with torsion, and to the case of nonminimal differ- 
ential operators. In this paper the method of ref. [6] is generalized to the case 
of operators of the type of the square root of the Laplace operator. In order 
that extraction of the root be meaningful, the operator under the sign of root 
should be nonnegative, i.e., eigenvalues should be only positive or zero. It 
is not an essential restriction as to the applicability of the method proposed 
because in physics we are mainly encountered with operators bounded from 
below. Up to our knoweledge, there is no available method in current litera- 
ture for the calculation of the DWSG coefficients for such type of operators. 
In Section 2 we compute the lowest E2 DWSG coefficient for the operator 
y^— + V{x) by using a method proposed with the expansion of the root. 
In Section 3 we generalize the method to the case of an arbitrary natural root, 
i.e., for the operator (V— -|- V{x)y^'^ where m is any natural and also to 
case of the operator \/—V^ + V{x) which cannot be represented as any power 
of the operator — + V{x). In Section 4, for the case of strictly positive 
operators under the sign of root, we propose a more simple method for the 
calculation of the DWSG coefficients. By using this method, we were able 
to shown that the problem of the calculation of the DWSG coefficients for 
strictly positive operators under the sign of root is exactly solvable. Namely, 
an explicit formula expressing the DWSG coeffcients for operators with root 
through those for operators without root is deduced. 
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2 Method for calculation of the DWSG coef- 
ficients for the operator of the square root 
of the Laplace operator 

We take as our space a compact n-dimensional Riemann manifold M with- 
out a boundary. The operators will act on the space of a vector bundle 
over the base M. The covariant derivative acting on objects with fiber (left 
understood) and base indices is defined by the rule 

k 
i=l 

where r^j^ and are the affine and bundle connections; = —^iuj'^^Tiab + 
iA^, and are the representation operators of local rotation group SO(n) 
under which 4'iJ.i...iJ,k transformed, uj^^ is a spin connection and is gauge 
potential. For the commutator of covariant derivatives we have 

k 

[V^, Vv\4>IJ.l...lJ.k ^ Rflil,u4'lJ.i...fMi-lXtMi+l...lJ.k + W^jUI/0jUl.../ifc, (2) 

i=l 

where W^,^ = 9^<^i/ — 81,00^ + [o;^, cUi,] is the bundle curvature, and the Rie- 
mann curvature tensor it!p^,^ is expressed through the affine connection F^^ 
as follows: 

pA Q -pA o -pA I -pA per pA per /o\ 

^pfiu ^fJ.^ pv '-'v^ pn ' (TH^ pv ^ av^ pp- K'-') 

As it is well khown [2,3,10], for a positive elliptic differential operator A of 
the order 2r there exists an asymptotic expansion of the diagonal matrix 
elements of the heat kernel exp(— tA) as i ^ 0+ in the following form: 

m 

where the summation is carried out over all non-negative integers m and 
Em{x\A) are the DWSG coefficients. 

In this section we consider a generalization of the method of [6] to the 
case of the operator 

A = ^-V2 + y(x), (5) 
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where V(x) is an arbitrary matrix with respect to bundle space indices. Fol- 
lowing the method [6], to obtain expansion (4) we use the representation of 
the operator exp{—tA) {t > 0) through the operator A resolvent 



c 

where the countour C goes counterclockwise around the spectrum of the 
operator A. For the matrix elements of the resolvent {A — X)~^ we employ 
the representation in the form 

<x\-^\x'>= [ S--'M^,x\k-\), (7) 

where l{x, x' , k) is a phase function and (t(,x, ,t', k/. A) is an amplitude [10, 11]. 
In the fiat space the phase l{x,x', k) = k^{x — x'Y is a Unear function of k 
and X for each x' . In the case of a curved manifold the real function 1 must 
be biscalar with respect to general-coordinate transformations, and must be 
a linear homogeneous function in k. The generalization of the linearity con- 
dition in X is the requirement for the mth symmetrized covariant derivative 
to vanish at the point x' with m > 2, i.e. 

{V^, V^, . . . V^^} 11=,. = [{v^, V^, . . . V^^} I] = 

/c^j for m = 1 and for m 7^ 1. (8) 

In eq. (8) the curly brackets denote symmetrization in all indices and the 
square brackets mean that the coincidence limit is taken. The local properties 
of the function 1 are sufficient to obtain the diagonal heat kernel expansion. 
The resolvent of the operator A satisfies the equation 

{A{x, V^) - A) r(x, x\ k- A) = -^5{x - x'), (9) 

yd 

and, therefore, in order to fulfill (9) it is sufficient to require that the ampli- 
tude cr(a;, x', k] A) satisfy the equation 

+ iV^O - A) (7(x, x\ k; A) = l{x, x'). (10) 
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The biscalar function If matrix with respect to bundle space indices 

and is defined by the conditions similar to eq. (8): 

[I] - 1, 

[{V^,V^,...V^^}I] = m>l, (11) 

the unity in eq. (11) is a matrix unity. To generate expansion (4), we 
introduce an auxiliary parameter e into eq. (10) according to the rule Z — > Z/e, 
A — > A/e, and expand the amplitude in a formal series in the powers of e 

oo 

a{x, x\ k;X)=Yl ^'^"''^m{x, x', k; A) (12) 

m=0 

(the parameter e then set equal to one). Then, eq. (10) gives us the recur- 
sion equations to determine the coefficients cr^, and, finally, this procedure 
leads to expansion (4) where the DWSG coefficients Ejn{x\A) are expressed 
through the integrals of [cr^] in the form [6]: 




Up to now we followed [6] very closely. Differences arise for the operator of 
the type of the square root of the Laplace operator when we are going to 
obtain the recursion relations for cr^. For the ordinary Laplace operator the 
recursion relations for follow directly from eq. (10) but it is not a case for 

the operator A = y^— + V{x) which we consider. Explicitly, the equation 
for a takes the form 

I °° 

(V V^^V^Z - ieVH - e^V^ - 2ieV^ZV^ + €^V{x) - A) ^ e"*(7„ = I (14) 

m=0 

We cannot just expand the square root of the operator in the powers of e in 
the Tailor series as in the case of the Laplace operator because V^/V*/ and 
the operator with e and do not commute and it is not clear in which order 
to place them in the Tailor formula. Therefore, to generate an expansion 
of the root in powers of e we first write down a general structure for the 
expansion of the root in the powers of e 

^Jv^lVH - ieVn - e2V2 - 2ieV^/V'^ + e^V{x) = ^V^IWH + e/i+ 
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e2/2 + . . . + + . . . (15) 

where fi, f2 and /,„ are to be found. In order to show how the method 
proposed works, first we compute the lowest Eq and E2 DWSG coefficients. 
To do this, we have to find the expansion of the root up to e^, i.e. we have to 
find only fi and /2. This can be done as follows: First, we take the square 
ofeq. (15) 

V^lVl - ieVH - e^V^ - 2«eV^/V^ + eV(a;) = V;,/V^/+ 



yJV^lVHeh + e/i yJV^lVH + e^fl + ^V^JV^e V2 + e V2 yJV^lVH + . . . ( 16) 

Then, comparing terms with the equal powers of e, we obtain the equations 
for /i and /2 



iVn - 2iV^lV^ = JV^lVHh + fJV^lVH, (17) 



- + V{x) = ft + ^V^lVHh + f2^JVf.lW^^l (18) 

From the left-hand side of eq. (17) it follows that a general structure of /i is 
/i = — ia^V — ib, where and b are ordinary vector and scalar functions, 
respectively, not operators. Substituting the general expression for /i into 
eq. (17), we get the following equations for and b: 



- 2iW^lW^ = -iJVJVna^V^' - ia^JVJVHV^ (19) 



- iyH = -i^V^,lVHh - ih^V^lVH - iai.VyVJV'l (20) 
From these equations we obtain 



V-/ V^/V'VV./V-/ 



(21) 



Finally, the equation for /2 takes the form 



- V' + V{x) = -(a^V^ + hf + yjv^lh + h^Jy^l (22) 
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Similarly to the case of /i we write down a general structure of /2 



/2 = Ci V' + 6-2^^^ + C^^u^^V" + C4 (23) 
Substituting (22) into eq.(24), we find 

1 



'2i?V2' 



(^2^ — 



2i? 2R i?V2 + 2i?V2 ^ 



V{x) a^V% I? V^R^/^ 1 Vf,R^'-V''R^/^ 



2i?V2 Ri/2 2i?V2 4i? 2i?V2^ 2i? 

2^ ^ SV2 ^ ^ 2R ' ^ ^ 

where R = V ^IVH. 

Thus, we have found the expansion of the root up to but it is obviously 
that it is possible in a similar way to find the expansion of the root up to 
any mth power of e because the equation for has a similar form to the 



equations for /i and /2, namely, \JV fjlVHfm + fmsj'^ ixl^^l- Consequently, 
writting down a general stucture of and defining /i, /2,..., fm-ii we can 
find the explicit expression for fj^ in the same way as it was done in the case 
of /i and /2. 

Having obtained the explicit expansion of the root 



yjv^l - ie(a^W + 6) + £^2 + ■ ■ ■ , (25) 
from eq.(14) we have the following equations for uo, (7i and (T2: 

(i?V2 _ A)ao = I, 



_ x)a, - i{ai_,V^ + b)ao = 0, 

(R^/^-X)a2-i{a^W+b)ai + (CiV^+C2^W+Cs^,V''V''+C4)ao = 0. (26) 
Prom eqn.(26) we obtain 



N = 



V{x) 



2A;2(VA;2 - A)3 2A;2(VA;2 - A)3 2Vk^{Vk^ - X)^ 



^ ^ (28) 

4(A;2)3/2(v/p _ A)3 4(A;2)3/2(v/p_ A)3' ^ ^ 

where we introduced the notation [V ■ ■ - ^ \i] = kalnu...\a (see [6]) and 
wrote down only terms which do not vanish after the substitution of the 
exphcit expression for l^i,,„x and the convolution with k'^k" . . .k^. 

Let us recall that the DWSG coefficients are given by (13) and we have 
to calculate the integrals in A and k. The integral in A is trivially calculated 
by using the residue theory, and, consequently, we obtain 



Eo{x) = J 



dJ'k 



fc2 



(27r)-y^ 



(29) 



(27r)-^ 



^ 4fc2 



4A;2 4(A;2)3/2 4(fc2)3/2 

To calculate the integral in k, we note that (see 

d^'k 



-)■ 



(30) 



1 Z*^ 
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= g^^k^K, (31) 

where fi'{^i/i2...Ai2s} ^^'^ symmetrized sum of metric tensor products. Inte- 
grating in A;^, we obtain 

2r(n) 

- (4.)V.r(n/2) ' (^2) 
1 , yjn) V{x)r{n - 1) l^:r{n) 



E2{X) 



(47r)"/2 ' 4r(n/2 + 1) r(n/2) 4r(n/2 + 1) 

)■ (33) 



- 1) llZnn - 1) 



4r(n/2 + l) 4r(n/2 + l) 
Using l^,^xa = - [6], we obtain 



^.W-7.JS^(?-V-W). (34) 



(47r) V2r(n/2 - 1) V 6 

Note that E2 obtained for the operator of the type of the square root of 
the Laplace operator coincides with E2 calculated for the Laplace operator 
up to a contant factor and, thus, the dependence of this coefficient on the 
space dimension is rather trivial (cf. with the case of nonmiminal operators 
[8] whose leading coefficients also are not a power of the Laplace operator). 
We will show in Section 4 that for strictly positive operators the same is 
true for the DWSG coefficients of an arbitrary order. It is an interesting 
problem whether it is also true for operators which have zero eigenmodes. 
The method proposed permits a generalization to the case of any rational 
root and can be also used for the calculation of the DWSG coefficients for 
the operator \/— + V{x) whose the square is not the Laplace operator. 

3 A generalization to the case of an arbitrary 
rational root and the operator which can- 
not be presented as a power of the Laplace 
operator 

In this section we generalize the method proposed to the case of an arbitrary 
rational root, i.e. for the operator (— + V(a;))^/'", where p and m are any 
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naturals. For the sake of simplicity, we actually consider the case of a natural 
root, i.e., the operator of the type (— + V{x)y^"^, where m is any natural 
number and show that the method can be easy generalized to the case of an 
arbitrary rational root. The equation for a has the form 

((-(V^ + ^V^/)(V^ + iV^'l) + - a) a{x, x\ k; A) = I{x, x'). (35) 

To generate the heat kernel expansion, we introduce an auxiliary parameter 
e into eq. (35) according to the rule / //e™''^, A ^ A/e and expand the 
amplitude in a formal series in the powers of e 

oo 

a{x, x\ k;\) = J2 e^+'^a,(a;, x', k; A). (36) 
We seek an expansion of the root in the form 

e"'V{x)y/"' = {V^lViy/"" + e"^/Vi + e"^/2 + . . . (37) 

Taking the mth power of eq. (37), we obtain the following equations for the 
unknown fi and 

m—1 

- iVH - 2iV^lV'' = Yl R'fiR"^'^'', (38) 

i=0 

m— 2 m— 2 m—1 

- + V{x) = R'fiR^^^fiR""'^'''^ + Yl R'f2R""~^~\ (39) 

i=0 j=i i=0 

where R = (V/i/V^/)^/"*. Note that in the case of an arbitrary ratioanl root, 
i.e., for the operators of the type (— + V{x)yl'^, we would be have the 
pth power of —i'V'^l — 22V^/V^ on the left-hand side of equation (37) that 
docs not present an untractablc problem for generating the expansion in the 
powers of e. All the same we again can write down a general structure in 
derivatives of the root and taking the mth power can find the unknowns in 
the expansion of the root. Thus, the method can be used in the case of an 
arbitrary rational root. Writting down general structures of fi and /2 

fi = -t{a^V^ + b), 
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f, = CiV^ + C2^V^ + Ca^.V'"'^^ + C4, (40) 

from eqn. (38) and (39) wc find the explicit expressions for /i and /2 and 
from eq. (35) obtain the recursion relations for cto, cti and (T2 

(R - X)ao = /, 
{R - A)(7i - i{a^V^ + b)ao = 0, 

{R - A)(72 - i(a^V'^ + b)ai + (CiV' + €2^^^+ 

Csm^V^V^ + C,)ao = 0, (41) 
where R — (Vj^ZVZ)^/"*. Taking the coincidence limits, we have 

[^0] = (A;2)l/m_A' 



^2(^2)2-2/m((^2)l/m _ ;)^)3 ^2(^2)2-2/m((^2)l/m _ ;)^)3 

V{x) 

m(F)i-V"»((A;2)iM_ A)2~ 
(m - l)k^kHZx (rn - l)ki'kH^i'^ 



^2(p)2-l/m((^2)l/m _ ^^3 ^2(p)2-l/m((p)l/m _ ;^)3" (^2) 

where we write down only terms which do not vanish after the substitution 
of the explicit expression for and the convolution with k^-^. 

Integrating in A and k, we obtain the lowest DWSG coefficients of the 
heat kernel expansion for the operator (— + V[x)Y/^ 

mVi^n) , , 



mr m(^ //? 
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Eq and E2 calculated in the case of an arbitrary natural root m coincide 
with Eq, E2 calculated in the particular case m=2 (see eqn. (32) and (34)). 
Note that in comparision with the case of nonmiminimal operators [9] the 
dependence on m and the dimension of space is rather trivial, namely, E2 
depends on m only through a constant factor. 

We now show how the method proposed works in the case of the calcula- 
tion of the DWSG coefficients for the operator \/— + V{x) which cannot 
be obviously represented as a power of the Laplace operator. We can use the 
old expansion (25) for the root ^V^IVH - ieVH - e^V^ - 2ieW^lVi^. Fur- 
ther, as usually, in order to generate the heat kernel expansion, we introduce 
an auxiliary parameter e according to the rule / — > //e, A ^ A/e and expand 
the amplitude cr in a formal series in powers of e 

00 

a,{x, x', k;\)=J2 ^'^""^m{x, x' , k; A). (45) 

m=0 

Then, the equations for ctq, (Ti, (J2 take the form 

{R'/' - X)ao = I, 



{R^'^ - A)(7i + {-i{a^V^ + h) + V{x))ao = 0, 



(i?^/' - A)(72 + {-iia^V^" + fe) + V{x))ai 

+(CiV2 + C2^,W + C^^.WV' + C4)ao = 0. (46) 



Prom eqn. (46) we find 

V{x) 



(V^- A)2' 



k,V^V{x) V\x) 

L^2Jne«, - _ ^^^2 + (^ _ A)3 ' ^^'^ 

where we write down only the terms with V{x)\ the terms with Z^ . q, coincide 
with that for the operator y^— -|- V{x). The first term in the expression 
for [(72] new vanishes after the integration in k because of an odd power of 
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k. Note also that in difference to the case of the operator y — + V{x) 

coefficient Ei is not equal to zero for the operator \/— V{x). 
Integrating in A and k, we obtain 

2T(n)V(x) 

= _ 2r(n)y^(x) 

Consequently, the entire E2 coefficient is 



Note that contrary to the case of the operator y — + V"(a;) the E2 coeffi- 
cient for the operator v— V{x) essentially depends on the dimension of 
space. We can also generalize the method proposed to the case of the opera- 
tor of the type (— V^)^/™ + V{x), where m is any natural number. Thus, we 
have shown that the method proposed can be modified and adopted for the 
calculation of the DWSG coefficients for various operators which involve the 
extraction of root and have calculated lowest E2 coefficient for three differ- 
ent operators. Amount of work needed to calculate the DWSG coefficients 
icreases very quickly with the growth of the order of the DWSG coeffcient in 
the method proposed. It is connected with rapid increase of the number of 
terms in the expansion of the root with the growth of the order of the DWSG 
coefficient. In fact, in the case where the operator under the sign of the root 
is strictly positive there exists a more simple and less laborious method for 
the calculation of the DWSG coefficients. 



4 More simple method for the calculation of 
the DWSG coefficients for strictly positive 
operators 

Let us again consider the operator A = \J + V{x). If this operator is 
strictly positive, i.e., it does not have any zero eigenmodes, we can used 
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instead of (6) the following representation: 



e-^^= ['^e-'^''\A-X)-\ (51) 
J 2tt 
c 

We demand that — + V{x) do not have any zero eigenmodes because in 
such a case the contour C can be drawn such that it encircles the whole 
spectrum of the operator —'V^ + V{x) and do not intersect anywhere the cut 
from infinity to zero along the negative half-axis which is needed in order that 
the extraction of root be meaningful. If the operator has zero eigenmodes, 
such a countour cannot be drawn because in this case we cannot draw the 
countour in such a way that it do not intersect the cut and simultaneously 
the contribution of eigenmodes be properly taken into account. By using this 
method, we can prove that the DWSG coefficients for operators with root are 
expressed through the DWSG coefficients for operators without root. Note 
that this method cannot be used for operators which have zero eigenmodes. 
However, the method with the expansion of root can be use also in the case 
of operators with zeromodes. Of course, for strictly positive operators both 
methods can be used and they yield coinciding results. 

Era coefficients in the method with the representation e~*^^^^ are given 
by the relation 

^-^"^'"^^ = / (4^ / '-^e-''^'[a^]ix,x,k;X), (52) 

where am{x, x, k; A) are the same as for the operator — + V{x). Using a2 
obtained in work [6] and calculating the integrals over A and k, we obtain 
the following E2 coefficient for the operator ^— + V{x): 

which coincides with E2 obtained in Section 2 by using the method with the 

expansion of root. 

Let us show by using the method with representation (51) that the 
DWSG coefficients for operators with root are explicitly expressed through 
the DWSG coefficients for operators without root. Let us consider the DWSG 
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coefficients for the operator — + V{x). According to [6], they are given by 
the relation 

E^{x\A) = I .An^ I '-^e-^[am]{x,x,k;X). (54) 



c 



Comparing it with (52), we see that the only difference is the power of A 
in the exponent, namely, it is equal to 1/2 in the case of the operator with 
root and 1 for the operator without root. We recall that a general term of 
[a{x,x,k; X)] has the following form: 

(F - A)« ' ^^^^ 

where F^^'^-^^' is expressed through the bundle curvature Wfj,^, and the Rie- 
mannian curvature tensor Rp^„. [cXmix , x , k; \)] is the sum of terms with 
various powers of a and s. It is very important for what follows that the 
difference a — s is fixed for the DWSG coefficient of a given order. This fact 
follows from the homogeneity property of the reccurent relations for (see 
[6]). a — s is equal to 1 + m/2 for the operator —'V'^ + V{x). Integrating over 
A and angles in n-dimensional space, we have for the DWSG coefficient in 
the case of the operator without root 



/ 



dtt-'+^'j,^^^^.,^-*- ^^e^"' (56) 



and for the operator with root 



/ 



where we have omitted common constant factors which coincide for two cases 
under consideration and have used formula (31). Integrating over k, we 
obtain r(^ + s — a + 2) — r(^^^) for the operator without root and 
2r(n — 2 + 2s — 2a + 4) = 2T{n — m) for the operator with root. It is 
very essential that the F-functions do not depend on a and s due to the 
homogeneity property and depend only on m. Therefore, the results obtained 
are true for any term in the expansion of cr^. Thus, the DWSG coeffients for 
the operator with root are expressed through the DWSG coefficents for the 
operator without root 

_ 2T{n -m) 
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where Emr are the DWSG coefficients for the operator with root. It is easy 
to check that Eq and E2 directly calculated in Section 2 for the operator 
Y^— + V{x) by using the method with the expansion of root (see formulas 
(31) and (34)) coincide with the DWSG coefficients given by the common 
formula (58). By using the representation with e^''^'' , similarly, it is easy to 
show that the DWSG coefficcnts for the operator with an arbitrary rational 
root, i.e., for operators of the type (— + V{x)y^'^, where p and q are 
any natural numbers, are expressed through the DWSG coefficents for the 
operator without root as follows: 

p q/pni/P^) p .... 

In the particular case of the square root, this formula yields (54) and in 
the case of natural root, i.e., p = 1, Eq and E2 given by (59) coincide with 
Eq and E2 explicitly calculated in Section 3, formulas (43) and (44). Thus, 
the problem of finding of the DWSG coefficients for operators of the type 
of rational root of a strictly positive operator is exactly solvable, i.e., the 
DWSG coefficients for operators with root are explicitly expressed through 
those for operators without root. 

Note that it would be of significant interest to calculate the DWSG coeffi- 
cients by using two methods proposed for an operator which has eigenmodes. 
Finding the difference between E^ obtained by two methods, we would be 
able to define the contribution of eigenmodes to the DWSG coefficients. 
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